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0\ . Abstract 

For the S-states of muonium and positronium, the hyperfine shifts to the order 

CJ ' a of a recently derived two-fermion equation with explicit CPT-invariance 

\my • are checked against the results of a nonrelativistic reduction, and the leading 

OO . a^ shifts are calculated. An additional hyperfine operator is discovered which 

can milden the singularity for r ^ of the Dirac hyperfine operator, such that 

the resulting extended operator can be used nonperturbatively. The binding 

^_ correction to magnetic moments is mentioned. 
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I. INTRODUCTION 



The energy levels of one-electron atoms are well described by the Dirac equation, par- 
ticularly when the nuclear hyperfine operators are included, and when radiative corrections 
/\ ■ are computed from the vector potential operator that is part of the electron momentum 

c^ ■ operator (Erickson 1977). But the Dirac equation applies also to muonium (e'/i"*"), where 

the "nucleus" is a structureless muon. A simple recoil correction reduces the electron mass 
mi by a factor m2/m = 1712/ {mi + ^2), where m2 is the muon mass. The resulting reduced 
mass /i„r = m,im2/m, is familiar from the two-body Schrodinger equation, but apart from the 
limit m2/m, ^ 1 — m,i/m,2 (Braun 1973), the correctness of this factor in the Dirac equation 
has not been proven. 

For positronium (e^e"*"), 1712/ m, = 1/2 excludes a perturbative treatment of recoil. In 
the past, precise calculations of muonium and positronium have used the Bethe-Salpeter 
equation, which treats both particles relativistically (Sapirstein and Yennie 1990). However, 
the complicated higher-order radiative corrections contain only few relativistic effects. An 
alternative method has been elaborated in which all relativistic effects are treated perturba- 
tively, by effective operators added to the nonrelativistic two-body Schrodinger Hamiltonian 
p'^/2finr + V{r) (Caswell and Lepage 1986). By this "nonrelativistic quantum electrodynam- 
ics method" (NRQED), the energy levels of muonium and positronium have been calculated 
to the order a^ (Pachucki 1997). To that order at least, all results of the Dirac equation with 
the above reduction factor are confirmed at the lowest nonvanishing power of /inr/""^, not 
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only in the coarse and fine structures, but also in the hyperfine structure. In the calculation 
of the muonium hyperfine splitting, the term a^^\^/m must also be included, but only in 
the approximation /x^^/m ~ m\/m2 (Kinoshita 1998). It is clearly more convenient to rely 
on the correctness of the Dirac expectation value of the relativistic hyperfine operator than 
to extract this term from four-photon exchange graphs in NRQED. 

Recently, a relativistic two-body equation with only eight components and explicit CVT- 
invariance has been constructed from the QED scattering matrix (Hackl et al 1998) for two 
structureless fermions (a short rederivation is presented in appendix A). Although both 
particles are treated symmetrically, the equation looks like a one-body Dirac equation, in 
which the mass mi and energy Ei of particle 1 are replaced by a relativistic reduced mass 
/i and reduced energy e: 

I2 = rmm2/E, e/i^={E^-ml-ml)/2mim2, (1) 

where E is the total cms energy (c = 1). This equation has predictive power for all values 
of Unr/'iTT', including the positronium value finr/fn = 1/4. For the nS-states, E is finr times 
a polynomial in /x^^/m, which has the order 2 to the order a^ and 3 to the order a^. In 
this paper, the second-order polynomial is checked against the NRQED results, and the 
third-order polynomial is constructed to an extent that may become relevant in the near 
future. 

For a simple understanding of the origin of the polynomial, the hyperfine operator is 
neglected in section 2. The energy levels are then obtained from those of the static Dirac 
equation, Ei/rrii = f{n*), by the replacement Ei/mi -^ e/ jj,, and by the subsequent evalu- 
ation of -E — m as a power series in / — 1. The resulting terms of order /i„r and fi'^j./m give 
directly the hyperfine-averaged energy levels, while the order ^\^/m? contains additional 
second-order hyperfine effects in the hyperfine-averaged levels, beginning at the order a^. 

The argument n* in f{n*) denotes the effective principal quantum number, n* = n + (5£, 
where 5(i is a (negative) quantum defect. With 5(i unspecified, the results of section 2 are 
more general than their derivation. This is so because the function 

f(n*) = (1 + aVn*')-^/2 ^ 1 - a^/2n*^ (2) 

solves not only the Dirac equation, but also the Klein-Gordon equation. The latter case has 
^^KG = [(^+1/2)^ — a^]^/^ — £ — 1/2, while the Dirac case has the orbital angular momentum 
i replaced by j = 1/2 for £ = and by j = ^ ± 1/2 for i > 0: 

siD=i-j-h i = [{j+ir-c^v'- (3) 

Thus expressions (1) and (2) occur also in the bound states of two spinless particles such as 
TT^K^ and tt^ti^. On the other hand, bound states of one lepton and one spinless particle 
(yU^7r+ or e~ ^He) require a more complicated equation with an asymmetric dependence 
on rrii and 7712, which produces the socalled Barker-Glover term at the order a^. (Atomic 
hydrogen requires also an asymmetric interaction, due to the proton's large anomalous mag- 
netic moment.) In section 3, the expectation values of the Dirac hyperfine operator of 
the two-fermion Dirac equation are expanded in powers of f{n*jy) — 1 to order a^fi^^/m^. 
The first-order expansion for Zeeman operators is mentioned, in which the CVT origin is 
particularly evident. 



Radiative corrections are ignored so far, but for the more complicated graphs, they may 
again be evaluated nonrelativistically. The simpler graphs can be evaluated by Dirac meth- 
ods also for the new equation. There are however other effective two-photon exchange op- 
erators that must be added in the two-fermion equation. For i > 0, they are 7 a^ / Qtt ^nrf^''^^ 
and Q;^L^/2yU,^^mr^. The first of these contributes at the order a^ /ti, where it is part of the 
socalled Salpeter shift. Its precise form for S-states is Ta^pQ/Gn^nr (Gupta et al 1989), 
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with 7 = e*" and C = Euler's constant 0.577. ..,with the prescription that expectation values 
of Pq are calculated by partial integration. This operator is mentioned here because its rel- 
ativistic analogue occurs in the hyperfine interaction at the order a^/i^^/m^, to be discussed 
in section 4. 

II. THE ENERGY LEVELS OF THE TWO-FERMION DIRAC-COULOMB 

EQUATION 

When the hyperfine operators are neglected, the resulting Dirac-Coulomb equation has 
exact solutions. In units of the relativistic reduced mass fi ([1|), the equation reads {h = c = 1) 
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with P = 7°, 75cri = ct = 7^7 as usual, and p^ = p/fi. 
The solutions are 
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with 7 defined in (^. With the definitions (|T]) of /i and e, the binding energy E — m is 
expanded in powers oi fr, — 1 as follows: 



E -m= ^Jm? + 2mim2{fD - 1) - m = finrifo - 1) 
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The precise expansion parameter is mim2{a / rfj^mY , but as n^ is also expanded in terms of 
a at a later stage, the mass dependence must be kept explicitly in (|^. To the order a^, one 
needs 
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It has been pointed out that the simple relation (|l|) between e/ ii and E'^ makes the expan- 
sion of E"^ — m? simpler than that oi E — m, but for the check against present NRQED it is 
necessary to expand E — m. In particular, one sees that for a^ the mass dependence ends at 
li^^/w?. This order in Unr/i^ receives contributions also from second-order hyperfine effects, 
which are discussed in the next section. 

In the frame of the two-fermion Dirac equation, all terms discussed so far originate from 
the one-photon exchange Born graph. Additional two-photon exchange graphs contribute 
already at the order a^/vr. These have only been evaluated with nonrelativistic approxima- 
tions, where most of them result in the operator 5{r), with matrix elements proportional to 
6eo/n^- The only two-photon exchange operator with a more complicated n-dependence is 
the "Gupta operator" (^), which has the following expectation value for £ = 
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with \1' = \I'(n) = dT{n)/dn. At this order in a, there is no finite-range two-photon exchange 
operator in the hyperfine structure. At the next order a^, there is no finite-range three- 
photon exchange operator for S-states in the Dirac-Coulomb equation, but there is one such 
operator to be added to the Dirac hyperfine operator, which is discussed in section 4. 
All terms in (|^) that require only the nonrelativistic values, /d — 1 = — a^/2ra^, (/^i — 1)^ = 
a^/4n^ etc. are valid for all combinations of spins and magnetic moments. At the order a^, 
there is only one such term, originating from the second term in the expansion (|^). It was 
originally discovered by Bechert and Meixner (1935). The complete series (^) was discovered 
for two spinless particles by Brezin et al (1970), with n*^ replaced by the appropriate Klein- 
Gordon value n*f^Q. It could not yet be confirmed experimentally. For positronium, on the 
other hand, the agreement between the NRQED calculation to the order a^ and experiment 
is not perfect, but it seems fair to say that all terms in the expansion (0) are now checked, 
for the extreme case finr/^^ = 1/4 (Pachucki and Karshenboim 1998, Czarnecki et al 1999). 

III. DIRAC HYPERFINE SPLITTING 

Inclusion of the hyperfine operator in the dimensionless reduced Dirac equation (^) leads 
to (Hackl et al 1998) 



f -V(p)-/?- 75(^1 + ^&^)P, 



^ = 0, ^« = -^o-i X o-2V^| (12) 



By standard perturbation theory, the first-order shift of e/fi caused by the hyperfine operator 
is (Rose 1961) 
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where f = j ± 1/2 is the total angular momentum, and K£, = 2{i — j){j + 1/2). Any small 
shift 5{e/ jj) = 6{E'^/2mim2) corresponds to a small shift 6E, 
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Thus the first-order hyperfine shift is 
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In the hyperfine structure of muonium, finr/in < oi makes the last term in ([15|) neghgible, 
while /d — 1 is required at most to the order ot^. Insertion of n^^^ ^ rr"^ + a^ /rr'{j + 1/2) 
gives 
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Turning now to flf£,, we find to order a^ and for j = £ + | 
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For j = \ (nS-states), the a^ correction in (^) becomes 
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For the "circular" states with n = j + ^ (which include the ground state) (H) gives n}) 



li fD = l/n, in which case fl/i^ is greatly simplified: 



(1) 



Jhf,D 



n — 



f-J 



wf + n/2 



f-J 



f + 1/2 n37(27 + 1)(27 - 1) / + 1/2 n^{2-f^ - 7) ' 



(20) 



This expression was used in the calculation of the muonium hyperfine splitting. For n = 1, 
the a^-component of the last bracket is 17q;'^/4, in agreement with (|TP|). It amounts to 
12 ppb in the muonium hyperfine splitting, Az/ = 4463 302 617(510) Hz (the error of 510 
Hz arises from the uncertainty of m^). As both C2(j = |) and its recoil correction (|15D 
are positive, inclusion of the latter one gives an (insignificant) increase of Az/. A recent 
experimental determination (Liu et al 1999) gives Au = 4463 302 776(51) Hz. 

The hyperfine splitting in P-states with total angular momentum / = 1 is complicated 
by the mixing of j = ^ and J = | states at the order a^; in this case the /i„r/'"^-expansion 
for a given power of a does not terminate (Pilkuhn 1995). Fortunately, the S-states are 
simpler. At the order a^, the n^Si states have a mixing between J = | and J = | (the 
S-D-mixing), and all S-states have contributions from the squares of the diagonal hyperfine 
matrix elements. Both effects are of second order in the hyperfine operator; their energy 
shifts {s/fi)hf may be calculated with Schrodinger wave functions and then simply added to 



the higher-order relativistic terms of {e/ jj)\^1. The S-D-mixing contributes —Aa^iJi\^/Qrn?n^ 
to the hyperfine splitting, which combines with +4a^/i3^/m2n^ from the second term of (pIBf ) 



into a total of 32a^/i^^/9m^n^, in agreement with the NRQED result (Pachucki 1997). This 
is presently the only confirmation of the expansion ([T5|) . The squares of the diagonal matrix 
elements may be calculated from an effective Schrodinger equation, which in units of fi is 
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where p comprises pz^ and all other interaction that may be approximated by a p -potential: 



Pl = -{dp + yp?+n^'+i)/p' 



The case at hand has 
To order a^, one obtains 
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where the term —3a'^6i'^finr/'^n'^ contains the desired Si\j-. Its contribution to the hyperfine 
splitting is — 16a^/i^^/3?7i^n^. Here, however, the NRQED results give only half that value, 
together with other terms, to be discussed in the next section. 

Before leaving the subject, we wish to propose a minor modification of the binding 
correction to the leptonic (^-factors that has been used in the measurement of the muon 
magnetic moment in muonium (Liu et al 1999). To order a^ and to second order in mi/m2, 
this correction has been calculated (Grotch and Hegstrom 1971, Faustov 1970) as 



9i^9i 



a 



nil 



1 ha' 

3 2m2 



2— 
m2 



(25) 



where Qi (i=l,2) are the ^f-factors of the two free leptons. The modification consists of 
replacing the mass factors by finr/'m- More precisely, we propose 
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The reason is that the Zeeman operator must be an odd function of E. In fact, any small 
change 6E caused by CPT-invariant perturbations must be odd in E. In the case at hand, 
CPT-invariance requires equal energy levels for muonium and antimuonium, even in the 
presence of a magnetic field. Both systems are described by a single eigenvalue equation, 
with E^ as eigenvalue (Malvetti and Pilkuhn 1994). First-order perturbation theory produces 



a small shift 5(-E^), from which 6E follows as in eq. ([l^ . For the hyperfine operator 



of equation ([T^), 6{E'^) goes as E~'^, due to jj,/E = mim2/E'^. For the Zeeman operators, 
S{E'^) is obviously independent of E, to order a^. For positronium, (^61) has been verified 
by a NRQED calculation (Grotch and Kashuba 1973). 



IV. THE GUPTA HYPERFINE OPERATOR 

According to the NRQED calculation to order a^, the total hyperfine splitting ^Ehj = 
£^(triplet)-i?(singlet) contains also the combination — 16a^/i^^(ln(a/ra) + \1' + C + 7/6 — 
l/2n)/3n^'m'^. The last term in the bracket may be decomposed into —l/n + l/2n, where the 
—1/n arises from the previously mentioned second-order hyperfine interaction. Comparison 
with (Q) shows that the n-dependence of the remaining part is identical with that of the 
"Gupta operator" . The complete a^-result for S-states is 

^eH^Inrqed = ^Efl^ + a4/!g + 8aVLi^;://3n3m^ (27) 

where AEU ^ contains all terms from the two-fermion Dirac equation, and 

A<,G = -f «^^(j>y2/.), F^, = | + ^ln^-21n2 + /.,, (28) 

where fhf is a numerical function of ^nr/^T^-, which is presently known in analytic form only 
for finr/'tTT' = and 1/4. The "Gupta" part of (pSD is the nonrelativistic expectation value 
of an operator that may be combined with the Dirac hyperfine operator in (|T^) : 

CT,f = a^^)[l + 2{^i/E)aHn{^p)]. (29) 

With (t\^J = —icTi X cT20il~i/Ep, the Gupta hyperfine operator goes as p^^ Inp, which is even 
more singular than the p~^ of the Dirac hyperfine operator. However, as it is only a first- 
order correction, it may be combined with the Dirac hyperfine operator into less singular 
forms, for example 

(Thf = —icTi X cr2Vmim2[E'^ — 2r7iim2a^ln(7p)]^^. (30) 

This combination goes as (plnp)~^ for p — > 0, which appears to admit a nonperturbative 
use in the Dirac equation. 

This extension is of little importance for QED bound states. But hyperfine operators 
appear also in quarkonium models, mainly in the form of Breit operators for heavy quarko- 
nium. For the vector (1~) and pseudoscalar (0~) mesons, one expects a small and constant 
hyperfine splitting in E'^, A = E'^{1^) — E^{0^) (Mannel 1998). However, a closer look 
reveils that A increases with decreasing meson masses, from 0.48 GeV^ for the heavy b 
quarkonium (Review of Particle Physics 1998) up to 0.57 GeV^ ~ ml for the p — it system, 
where NRQED expansions would diverge. Moreover, m^ Jm^ = 30 excludes a perturbative 
treatment of the hyperfine operator. It is true that the quark model for pions must differ 
drastically from any QED analogue, for example in its dependence on mi and m2. But 
the pion is only the lightest meson in a long list of light mesons whose quantum numbers 
all agree with the naive quark model. Until the masses of these mesons are calculated by 
fundamental methods such as lattice QCD, it may be allowed to replace the nonrelativistic 
QCD Breit operators by the relativistic QCD hyperfine operator. This operator has a similar 
structure as in QED. Even in the absence of a detailed quarkonium model, the experimental 
increase of A hints at an increase of the hyperfine operator for small E"^, in agreement with 
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APPENDIX A: REDERIVATION OF THE TWO-FERMION EQUATION 

The standard 16-component formalism for two fermions contains two kinetic energy 
operators otiPi (i=l,2), with a.i = jfcTi and (Ti = Pauli spin matrices. The essential point 
of the eight-component formalism in the cms {pi = — Pg = p) is the absence of C!t2P2 = 
— 7|cr2p in the free two-body equation, to which the interaction operator is added. The 
desired equation is easily constructed from the kinematical constraints: The direct product 
V'o = ''/'i ® "^2 of two free particle solutions tpi and 1^2 satisfies two Klein-Gordon equations, 
which are required in the cms, where one may use pi = p\ = p^: 

{pf - m\),p, = {pf - rriDij, = p^^Jo- (Al) 

After elimination of Pi— P2 = idn — idt2 using (p^ — P2) V'o = ij^i — ""^2) (Pi + vVl'^'^Oi there 
remains a single equation containing p\+ p\^ E: 

(A;2 - p2)V'o = 0, AE'^k'^ = E^ -2E'^{m\ + ml) + {ml-mlf. (A2) 

The equation applies to particles of any spins. For fermions, it must be linearized in p. Here 
it suffices to write p^ = (tsctip)^, with 7^ = 1 and one set of Pauli matrices, cti. At the 
same time, k'^ must be factorized, writing 

k^ = {e-fif3){e + fi/3), /3^ = 1. (A3) 

For two structureless fermions, the symmetry of k"^ under the exchange ml ^^ m\ must be 
preserved in the factorization, which leads to the expressions ([^) for /x and ej 11. In a last 
step, the E is removed from the denominators of yU and e: 

{Ee- EfiP -'j5criEp)tjJo = 0, Ee = -{E^ - ml - ml), Efx = mim2. (A4) 

The factor E in front of j? = — iV is absorbed by a rescaling of the variable r, after which 
( |A^ ) is an explicit eigenvalue equation for E'^. A slightly more convenient dimensionless 



In the 16-component formalism, the direct product ipQ = ipiD ® ^^20 of two free Dirac 



form is obtained by dividing by Efi = mim2 as in equation (^. 

In the 16-component formalism, the ( 
spinors satisfies two free Dirac equations, 

(pO-7,V.P.-m,A)4''^=0, A = 7°, Pl = ^du. (A5) 

The sum of these equations in the cms gives, with p^ = —p2 = p, 

[E - (7iVi - 72^2)^ - m^|3^ - m2/32]V'f ^ = 0. (A6) 



In the following, ( [Aq ) is transformed into ( |A^ ) before the corresponding interaction oper- 
ator is added. The motivation for this step comes from the details of the interaction, but 
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the main point is easily stated: In a first approximation, the differential approach with 
interaction included replaces the E in (A6) hy E — V^^\ where V^^"^ = —a/r is the main 



part of the Fourier transform of the first Born approximation. When this form is reduced 
to (|A4|) , then E'^ in Ee is replaced by (E — V^^^y. The second-order Born approxima- 
tion provides several additional operators, which in leading order cancel the squares of the 
first-order operators, such as the {V^^^y in {E — V^^^y. These cancellations occur not only 
in the differential equation scetched here, but also in the Bethe-Salpeter equation and in 
NRQED calculations. When the interaction is added in the eight-component form ([A4|) , E"^ 



is automatically replaced by E"^ — 2EV^^^; the operator (V^^^y is absent. In this sense, the 
first Born approximation in the eight-component scattering includes the leading terms of 
the second Born approximation in the 16-component scattering. 
To achieve the reduction from ( |A6| ) to (^), ipo is divided into two octets i/jqp and 



Xop, which have 7^ = 7! = 75 and 7^ = — 7I = 75, respectively. The round bracket is 
75(ci ~ (^2) = 75^cr when acting on ipQp, and 75(cri -|- (T2) = ^cr when acting on Xop- In 
the chiral basis, 7^ and 7I are diagonal: 

5_/-LU\ a _ I ^ ^ \ I _ I i'ir] I _ ( i'rr 



where the indices r and / (= righthanded, lefthanded) refer to the eigenvalues ±1 of jf. 
Each Pi exchanges ir with il; P2 exchanges ipop with xop, while /3 = /5i/32 exchanges each 
upper quartet with the lower one. 

^2^op = Xop , P2X0P = ^op h^ + Ts/? = 0. (A8) 

Consequently, ( |A6| ) is decomposed as follows: 

(E - 75pAcr)?/;op = "^+Xop, (^ - 75Pcr)Xop = "^+^0P, m± = m2±(3mi. (A9) 

Using the first equation for the elimination of Xopi 0116 obtains for the second one 

{E - -f5pa){m+)-\E - 75pAcr)^op = m+ipop. (AlO) 

Multiplying this equation by m_|_ and using 

m+75 = 75m_, (pcr)(pAcr) = 0, (All) 

one arrives at the following equation 

[E^ - Ej5{pcrm_/m+ + pAa) - m^^op = 0. (A12) 

As a rule, an elimination of components produces operators of second order in p. This is 
prevented here by ( |A11 ). The factor m^/m^ is removed by the transformation 



V'op = cipo, c "^75 = 75C, cere = crm^/m_ , cAcrc = A(t. (A13) 

c = (m+m_)"^/2 [m2 + mi/5(l + ctiCT2)/2] = (m+m_)~^/2(m+ - 2miA,/3), (A14) 



where A^ = (1 — cricr2)/4 is the projector on singlet spin states. The inner bracket in ( |A12|) 



becomes 2pcr^, and division by 2 produces (|A^). As mentioned before, the result is trivial, 
but the c-transformation will be needed again for the S-matrix, S* = 1 + iT . The 16 x 16 
form of the fermion-fermion T- matrix in the cms is Tij = u[u2TuiU2, where Ui and u[ are 
the free Dirac spinors for the in- and outgoing fermions, ipi = Uie'^'^' , (pi = kiVi — Eit. In the 
{ip, x) basis, analogous expressions are defined for the ingoing ipop and xop, 

iJop = ve"^, XoP = we"^, = 0i + 02, (A15) 

and for the outgoing ip'^p, Xop- '^if is now expressed as f'^T^f, f'^T^^w, w'^^T^^^v and 
w'^T^w. The first Born approximation T>j has T^„, = T^^, = 0. The product of Dirac 
matrices 7f72/x is /?i/32(l — 11(^112^2) , with 7j'7| = +1 in the t>-components and —1 in the 
w-components, respectively. Remembering u[ = u'^ (3i and /?^ = /3| = 1, one finds 

T^f = -ATiaq-^[v'\l - ctiCT2)v + w'\l + (Ti(T2)wI (A16) 

with q = k — k' [q^ = Q), where k and k' are the in- and outgoing momenta of particle 1. 
In general, T^w and Tyj^ appear only for an odd number of matrices jSi. 

One-loop graphs contain two fermion propagators, the product of which may be written 
symbolically as 

P = (^ - m)"^X(/ - m)-^ = {i + m)X(^' + m')/(/ - m^){p'^ - m'^), (A17) 

where X may be any operator. When both propagators occur on one fermion line i as in 
radiative corrections, one has m = m' = rrii otherwise m = mi, m' = m2. The product 
( |A17| ) may be decomposed as follows: 

P = P++-P__, P++ = 2(^X/ + mm'X)/(p2-m2)(p'2-m'2), (A18) 

P = (|( - m)X(j^' - m')/{p'^ - m^){p'^ - m'^) = (^ + m)"^X(/ + m')-^ (A19) 

For bound states, P is of higher order in a because it contains a product of two antifermion 

poles. The leading radiative corrections contain no /3, just as ( |A16|) . Two-photon exchange 

is linear in /5, because ;^Xf contains /3i/92 = /? in this case. 

Tif can always be put into the form w'^^Mv by means of ([X9|). For P ^ 0, one finds 

M = Tu,m~^{E - 75fcAcr) + {E - -f5k'(T)m-^T^. (A20) 

Before insertion into ( [A4| ) , M must also be c-transformed. 

The form w'^^Mv is no longer hermitian. It may be compared with the 2x2 single-fermion 
scattering matrix T^f = u'l MsUr, which is also complete and nonhermitian. 
The first Born approximation ( |A16|) gives 

m+M^^) = -ATTaq'^2{E - i-f^k cti x ^2). (A21) 

Notice the absence of k' . Its Fourier transform produces the interaction operators, which 
in the dimensionless variable p = fir lead to equation ([l2|) . Thus the only addition to 
the almost trivial Dirac-Coulomb operator (|^) is a hyperfine operator, which is not totally 
unexpected either. Apart from the replacement m^^ -^ m^^, which was already found by 
Breit, it differs from the standard hyperfine operator in two respects: the hermitization has 
been "forgotten", and the dimensionless form has ?7i~^ replaced by E~^. 
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